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Abstract.  The  stable  matching  problem  has  many  applications  to  real 
world  markets  and  efficient  centralized  algorithms  are  known.  However, 
little  is  known  about  the  decentralized  case.  Several  natural  randomized 
algorithmic  models  for  this  setting  have  been  proposed  but  they  have 
worst  case  exponential  time  in  expectation.  We  present  a  novel  structure 
associated  with  a  stable  matching  on  a  matching  market.  Using  this 
structure,  we  are  able  to  provide  a  finer  analysis  of  the  complexity  of  a 
subclass  of  decentralized  matching  markets. 
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1  Introduction 

The  stable  matching  problem  and  its  variants  have  been  widely  studied  due  to 
real  world  market  applications,  such  as  assigning  residents  to  hospitals,  women 
to  sororities,  and  students  to  public  schools  [1-3].  In  a  seminal  paper.  Gale  and 
Shapley  hrst  proposed  an  algorithm  to  find  a  stable  matching  in  the  basic  two- 
sided  (bipartite)  version  [4] .  Others  have  subsequently  investigated  the  structure 
of  the  set  of  stable  matchings  [1] .  However,  most  prior  work  involves  centralized 
algorithms  to  find  stable  matchings  where  the  entire  set  of  preferences  is  known 
to  some  central  authority.  In  some  cases  the  algorithms  are  not  totally  central¬ 
ized,  but  the  participants  are  subject  to  strict  protocols  where  only  one  side 
of  the  market  can  make  proposals.  Nevertheless,  many  applications  of  stable 
matchings  have  no  central  authority  or  enforcement  of  protocols,  such  as  college 
admissions  and  the  computer  scientist  job  market.  Therefore  we  investigate  this 
problem  in  a  decentralized  setting,  where  members  of  both  sides  of  the  market 
can  make  proposals. 

One  major  open  question  in  decentralized  stable  matching  concerns  whether 
natural  and  efficient  algorithms  exist.  To  this  end,  Yariv  argues  that  natural 
processes  will  find  a  stable  matching  and  provides  experimental  support  [5]. 
Roth  and  Vande  Vate  propose  a  class  of  randomized  algorithms  to  model  the 
decentralized  setting  and  show  that  algorithms  in  this  class  converge  to  a  stable 
matching  with  probability  one  [6].  At  each  step  these  algorithms  match  two 
participants  who  form  a  blocking  pair  (who  prefer  to  be  matched  with  each  other 
over  its  partner)  of  the  current  matching.  However,  they  present  no  expected  time 
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complexity.  Ackerman  et  al.  investigate  one  particular  algorithm  in  this  class, 
the  better  response  algorithm  (or  random  better  response  dynamics).  In  each 
step  of  this  algorithm,  one  blocking  pair  is  chosen  uniformly  at  random.  For  this 
algorithm,  they  show  worst  case  instances  that  take  exponential  time  to  reach  a 
stable  matching  in  expectation  [7]. 

Since  the  better  response  algorithm  is  natural  but  takes  exponential  time  in 
the  worst  case,  can  we  find  a  natural  subclass  of  matching  markets  which  do 
not  require  exponential  time?  Ackermann  et  al.  show  that  the  better  response 
algorithm  only  requires  polynomial  time  for  one  class  of  problem  instances,  those 
with  correlated  preferences  [7].  However,  correlated  preferences  require  that  a 
participant  obtains  the  same  benefit  from  a  partnership  as  its  partner.  This 
significantly  limits  the  preference  structures  allowed  in  the  matching  market. 
Therefore,  we  investigate  other  structural  properties  of  stable  matching  markets 
which  facilitate  faster  convergence. 

In  this  paper  we  make  progress  toward  answering  the  previous  question  by 
expanding  the  subclass  of  markets  with  polynomial  time  convergence  guaran¬ 
tees.  For  this  purpose  we  associate  a  directed  graph,  called  the  jealousy  graph, 
with  each  stable  matching.  It  turns  out  that  this  structure  is  a  key  factor  in  de¬ 
termining  the  convergence  time  of  the  better  response  algorithm.  The  jealousy 
graph  is  a  directed  graph  where  a  vertex  v  corresponds  to  a  pair  in  the  stable 
matching  and  an  edge  {u,  v)  is  present  if  one  member  of  the  pair  v  prefers  a 
member  of  the  pair  u  to  its  partner  in  the  stable  matching.  The  strongly  con¬ 
nected  component  graph  of  this  jealousy  graph  provides  a  decomposition  for  that 
stable  matching.  Our  intent  is  to  formalize  a  notion  of  structure  using  jealousy 
graphs  and  the  corresponding  decompositions.  In  particular,  we  find  that  the 
strongly  connected  components  of  this  graph  give  insight  into  the  complexity  of 
that  market.  Gusfield  and  Irving  provide  a  structural  property  of  stable  match¬ 
ings  which  describes  the  set  of  stable  matchings  and  the  relation  between  them, 
whereas  our  structures  relate  to  individual  stable  matchings  and  the  distributed 
process  by  which  these  stable  matchings  are  achieved  [1]. 

With  a  decomposition,  we  associate  a  size  and  depth.  Our  main  result,  The¬ 
orem  24,  states  that  for  a  matching  market  of  size  n  with  a  decomposition  of  size 
c  and  depth  d,  the  convergence  time  is  Therefore,  for  constant 

size  and  depth  decompositions,  we  demonstrate  that  the  better  response  algo¬ 
rithm  requires  only  polynomial  time  in  expectation  to  converge  for  an  expanded 
class  of  matching  markets.  This  indicates  that  the  jealousy  graph  and  decom¬ 
position  structures  partially  answer  the  convergence  questions  of  the  decentral¬ 
ized  stable  matching  problem.  As  an  application  of  our  work,  we  demonstrate 
how  Theorem  24  provides  theoretical  justification  for  the  simulated  results  of 
Boudreau  [8].  We  also  conjecture  that  these  structures  will  provide  a  means  of 
predicting  which  stable  matchings  are  likely  to  be  achieved  when  there  are  mul¬ 
tiple  stable  matchings,  a  question  that  others  in  the  literature  have  investigated 
[9,5,10]. 

In  the  remainder  of  this  section  we  formalize  our  model  and  present  the 
basic  concepts.  In  section  2,  we  present  some  useful  structural  properties  of 
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the  jealousy  graph  and  decomposition.  In  section  3,  we  have  our  convergence 
result.  Section  4  contains  the  application  of  our  work  to  [8],  and  section  5  is  our 
conclusion. 

1.1  Basic  Stable  Matching  Concepts 

We  start  with  the  basic  definitions  of  matching  markets  and  stable  matchings. 
Those  familiar  with  the  matching  literature  will  notice  that  we  restrict  prefer¬ 
ences  to  be  complete  and  strict. 

Definition  1  (S',  P)  is  a  matching  market  if  S  =  for  some  disjoint  sets 

M,  W,  \M\  =  \W\  and  P  =  where,  for  s  €  M,  is  a  total  order  over 

W,  and  for  s  €  W,  yg  is  a  total  order  over  M. 

We  say  a  matching  market  has  size  n  if  |M|  =  |IT|  =  n. 

Definition  2  A  matching  on  the  set  S  is  a  function  fx  :  S  ^  S  such  that'is  €  S, 
ix{ix{s))  =  s,  s  G  M  ^  ij,{s)  €  ITlJ{s}  and  s  GW  ^  ix{s)  G  MlJ{s}. 

We  say  that  a  participant  s  G  S  is  unmatched  by  a  matching  /r  if  /i(s)  =  s. 
We  also  assume  that  all  participants  prefer  to  be  matched  to  anyone  than  to  be 
unmatched.  Observe  that  fj,  can  be  thought  of  as  a  collection  of  pairs  (m,  w)  if 
we  allow  self  loops  (s,  s)  for  unmatched  participants. 

Definition  3  A  matching  on  the  set  S  is  a  perfect  matching  if  fj,{s)  ^  s  for  all 
s  G  S. 

Given  a  matching,  if  there  were  a  man  and  a  woman  who  each  preferred  the 
other  to  their  partner,  this  causes  the  matching  to  be  unstable.  Therefore  any 
stable  matching  must  have  no  such  pairs.  We  call  such  a  pair  a  blocking  pair, 
defined  formally  here: 

Definition  4  Let  (S,  P)  be  a  matching  market  and  p  be  any  matching  on  S. 
A  blocking  pair  for  p  in  (S,  P)  is  a  pair  (jn,  w)  such  that  m  G  M ,  w  G  W, 
p{m)  ^  w,  w  ym  m(w),  and  m  y^  p(w). 

Definition  5  Let  (S,  P)  be  a  matching  market.  A  matching  p  on  S  is  a  stable 
matching  for  (S,  P)  if  it  has  no  blocking  pairs  in  (S,  P) . 

The  following  three  concepts  will  be  useful  since  we  will  deal  with  subsets  of 
the  matching  market. 

Definition  6  A  balanced  subset  of  a  matching  market  {S,P),  S  =  M^jW,  is 
a  subset  S'  C  S  such  that  |S''p|M|  =  |S''p|IT|. 

Definition  7  A  matching  p  is  locally  perfect  on  a  balanced  subset  S'  C  S  if 
h'iS')  C  [S')  and  pIs’  is  o.  perfect  matching  on  S' . 

Definition  8  Let  p  be  a  stable  matching  on  a  matching  market  [S,  P).  A  match¬ 
ing  p'  is  p-stable  on  a  balanced  subset  S'  C  S  if  p'  ]s'=  h']  S'- 
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1.2  Better  Response  Algorithm 

The  class  of  algorithms  introduced  by  Roth  and  Vande  Vate  [6]  involve  randomly 
choosing  a  blocking  pair  of  the  current  matching  and  creating  a  new  matching 
by  matching  the  participants  in  the  blocking  pair  with  each  other.  This  resolves 
the  chosen  blocking  pair. 

Definition  9  A  blocking  pair  (x,y)  in  a  matching  fj,  is  resolved  by  forming  a 
new  matching  y!  where  y' (x)  =  y,  y' (y(x))  =  y(x)  if  y{x)  x,  y! {y{y))  =  y{y) 
if  h{y)  y,  and  y'{s)  =  y{s)  for  s  ^  {x ,  y ,  y{x) ,  y{y)} . 

This  process  is  repeated  until  a  stable  matching  is  reached.  The  better  re¬ 
sponse  algorithm  defined  in  [7],  is  the  algorithm  in  this  class  where  the  blocking 
pair  is  chosen  uniformly  at  random  from  all  blocking  pairs  of  the  current  match¬ 
ing.  Note  that  this  algorithm  results  in  a  sequence  of  matchings.  A  valid  sequence 
of  matchings  is  any  sequence  where  each  matching  is  formed  by  resolving  one 
blocking  pair  in  the  previous  matching. 

We  focus  on  the  better  response  algorithm  since  the  uniform  distribution 
on  blocking  pairs  facilitates  our  analysis  and  we  believe  it  provides  insight  into 
the  more  general  class  of  algorithms.  This  algorithm  also  serves  as  a  model  of  a 
distributed  stable  matching  market. 

1.3  Jealousy  Graph  and  Related  Definitions 

In  order  to  analyze  matching  markets,  we  represent  the  preference  structure  as 
a  directed  graph.  While  we  lose  some  of  the  preference  information,  we  retain 
critical  relationships  relative  to  the  stable  partners.  In  section  3  we  will  provide 
bounds  on  convergence  based  on  this  simpler  structure. 

Definition  10  The  jealousy  graph  of  a  stable  matching  y  on  a  matching  market 
{S,P)  is  defined  as  the  graph  =  {V,E)  where,  for  each  pair  {x,y{x)},  x  G 
S,  there  is  a  vertex  W{a;.^(a:)}  G  V  and  E  =  {(u{x,y}, "ytai'.y'}  G 
V,  and  either  x  )^yi  x'  or  y  y'}. 

The  jealousy  graph  can  provide  insight  into  the  complexity  of  stabilization.  For 
example,  suppose  the  jealousy  graph  for  a  stable  matching  y  is  one  large  clique. 
Even  when  all  but  one  pair  of  the  participants  are  matched  with  their  partner  in 
y,  there  are  still  many  blocking  pairs.  Therefore,  the  better  response  algorithm 
would  be  unlikely  to  choose  the  blocking  pair  that  would  result  in  a  stable 
matching.  This  greatly  hinders  convergence  to  the  stable  matching. 

On  the  other  hand,  suppose  the  jealousy  graph  for  ^  is  a  DAG.  Then  there 
is  at  least  one  vertex  with  no  incoming  edges.  This  means  each  partner  in  the 
corresponding  pair  is  the  other’s  first  preference.  Consequently,  this  will  remain 
a  blocking  pair  until  it  is  resolved,  so  we  would  expect  such  a  pair  to  be  resolved 
in  O(n^)  time  under  the  better  response  dynamics.  Moreover,  once  resolved,  the 
match  will  remain  unbroken  since  neither  partner  will  ever  be  involved  in  any 
blocking  pairs.  Ignoring  this  pair  will  result  in  at  least  one  other  source  vertex 
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of  the  graph.  Inductively,  these  pairs  will  be  resolved  in  O(n^)  expected  time. 
This  results  in  an  expected  convergence  time  of  0{n^)  for  the  matching  market. 
It  should  be  noted  that  the  class  of  correlated  markets,  for  which  Ackermann  et 
al.  prove  the  better  response  algorithm  requires  only  polynomial  time,  falls  into 
this  special  case. 

When  it  is  a  DAG,  the  jealousy  graph  provides  an  order  in  which  the  pairs 
will  likely  be  resolved  to  reach  /x,  namely,  a  topological  sorted  order.  However,  a 
matching  market  might  not  fall  into  this  extreme  case  as  there  could  be  cycles 
in  the  jealousy  graph.  Therefore,  we  define  a  decomposition  which  is  a  DAG 
obtained  from  the  jealousy  graph. 

Definition  11  Let  be  the  jealousy  graph  of  a  stable  matching  fj,  for  a  match¬ 
ing  market  {S,P).  A  fj,- decomposition,  is  a  graph  of  components  of  such 
that  if  u,  V  are  in  the  same  strongly  eonneeted  component  of  then  they  are  in 
the  same  eomponent  in  Pf,,  and  if  edge  {A,  B)  is  in  then  there  is  a  path  from 
a  vertex  in  A  to  a  vertex  in  B  in  J^. 

We  call  the  strongly  connected  components  of  stable  components.  Observe 
that  is  a  directed  acyclic  graph.  Therefore  it  induces  a  partial  order  on  the 
stable  components.  Sometimes  it  will  be  simpler  to  refer  to  the  decomposition 
as  Pfi  =  (17,  where  7T  is  a  partition  of  S  into  sets  corresponding  to  the  stable 
components  of  p^  and  ^  is  the  induced  partial  order  on  those  components.  As 
a  slight  abuse  of  notation,  we  will  use  the  term  stable  component  to  refer  to 
both  the  connected  component  in  the  decomposition  and  the  set  of  participants 
corresponding  to  this  component. 

In  dealing  with  partial  orders  we  will  use  the  concept  of  a  downset.  A  downset 
of  a  partially  ordered  set  77  with  partial  order  ^  is  any  set  such  that  for  A,  77  G  77, 
if  A  is  in  the  set  and  B  <  A,  then  B  is  in  the  set.  The  downset  of  an  element 
A  G  77  is  Down{A)  =  {B\B  ^  A}.  When  the  elements  of  77  are  sets  themselves, 
as  in  the  case  of  decompositions,  we  will  denote  union  of  sets  in  Down{A)  as 
D(A)  =  UBer>ou)n(A) 

For  our  complexity  results  we  need  the  following  two  notions: 

Definition  12  The  depth  of  a  stable  component  A  of  a  ^-decomposition,  p^,  is 
the  length  of  the  longest  path  in  p^  from  any  souree  vertex  to  va  ■  The  depth  of 
p^  is  defined  as  max^gp^  depth(A). 

We  will  say  that  a  stable  component  A  is  on  level  j  if  depth{A)  =  j.  Minimal 
stable  components  are  on  level  0.  Intuitively,  we  would  expect  components  on 
lower  levels  to  converge  to  the  stable  matching  sooner  than  those  on  higher 
levels. 

Definition  13  The  size  of  a  p-decomposition,  Pp,  is  defined  as  maxAep,^  size{A). 

Intuitively,  components  with  smaller  sizes  can  have  less  internal  thrashing  so 
they  will  converge  to  the  stable  matching  more  quickly  than  larger  components. 
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2  Structural  Results 

2.1  Any  Digraph  can  be  a  Jealousy  Graph 

These  structural  notions  would  not  be  very  enlightening  if  all  matching  markets 
had  the  similar  jealousy  graphs  and  decompositions.  However,  the  following  re¬ 
sult  shows  that  any  directed  graph  is  the  jealousy  graph  associated  with  a  stable 
matching  for  some  matching  market.  The  proof  can  be  found  in  the  Appendix. 

Theorem  14  Given  any  directed  graph  G  with  n  vertices,  there  is  a  set  S  = 
{mi,Wi\i  =  1,2, ...  ,n}  and  preferences  P  =  {Pmii'^wi  |1  <  *  <  fi}  such  that 
(S',  P)  is  a  matching  market  with  a  stable  matching  fj,  where  p.{mi)  =  Wi  and 
J>.  =  G. 

2.2  Properties  of  Decompositions 

In  this  section  we  prove  several  structural  properties  of  the  jealousy  graphs  and 
decompositions  essential  to  our  main  convergence  result.  The  proofs  are  in  the 
Appendix.  The  first  property  says  that  if  there  is  a  path  from  one  vertex  to 
another  in  the  jealousy  graph,  then  the  first  vertex  must  be  in  the  downset  of 
any  component  containing  the  second  vertex. 

Lemma  15  Given  a  matching  market  (S,  P)  with  a  stable  matching  fi,  let 
be  the  jealousy  graph  associated  with  fj,.  Let  v^rn.w}  ^iT-d  be  vertices  in 

J^.  Suppose  v^rn' ,w'}  G  ^  for  a  stable  component  A  of  a  fj,- decomposition  = 

(7T,  ^).  If  there  is  a  path  from  to  V[rn',w'}>  then  m,w  €  D(A). 

Using  this  lemma,  we  prove  that  no  member  of  a  stable  component  can  prefer 
anyone  outside  of  the  downset  of  that  component  to  his  stable  partner. 

Lemma  16  Given  a  matching  market  {S,  P)  with  a  stable  matching  p,,  let  = 
[n,  he  a  p-decomposition.  For  AgII,  a€  A,  s€S  —  D(A),  p{a)  )^a  s. 

A  further  property  is  that  if  there  are  two  stable  matchings  with  distinct 
decompositions,  the  intersection  of  the  downsets  of  stable  components  must  be 
mapped  to  itself  in  both  stable  matchings. 

Lemma  17  Given  a  matching  market  {S,P)  with  stable  matchings  p,p' ,  let  p^ 
and  Pfj_>  be  respective  decompositions.  Let  A  6e  D(A)  for  some  stable  component 
X  of  Pfj,  and  B  he  T>Y  for  some  stable  component  Y  o/p^/.  Then  p{A(^B)  = 
p'{A(}B)=A(}B. 

Our  hnal  result  shows  that  forming  a  stable  matching  on  the  downset  of 
a  stable  component  cannot  increase  the  size  or  depth  of  the  decomposition  of 
another  stable  matching. 

Lemma  18  Given  a  matching  market  {S,P)  with  stable  matchings  p,  p' ,  let  p^ 
and  Pfj,'  be  respective  decompositions.  Suppose  the  size  of  p^  is  c  and  the  depth 
is  d.  Let  A  be  a  stable  component  of  p^i .  Then  there  is  a  stable  matching  p" 
such  that  p"  ]r)^,{A)  =  h''  ]'d^,(A)  ond  p"  1s-d,^/(A)=  M  1s-d^/(A)-  There  is  also 
a  p" -decomposition  on  S  —  D^/(A)  of  size  at  most  c  and  depth  at  most  d. 
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3  Convergence 

In  this  section  we  prove  our  convergence  result.  The  proof  uses  two  main  ideas. 
First,  in  the  following  sequence  of  lemmas,  we  show  that  a  stable  component 
will  converge  to  a  locally  perfect  matching  in  time  that  is  only  polynomially 
dependent  on  the  size  of  the  entire  market.  Then  the  proof  of  Theorem  24  uses 
this  to  bound  the  time  it  takes  for  all  components  of  the  decomposition  to  reach 
a  stable  matching. 

For  this  section  we  will  assume  (S',  P)  is  a  matching  market  of  size  n,  /i  is  a 
stable  matching  on  S,  and  (7T,  be  a  /r-decomposition. 

The  following  lemma  says  that  if  a  matching  is  not  locally  perfect  on  a  stable 
component  of  a  /i-decomposition,  then  there  is  a  blocking  pair  which  is  in  fj, 
between  two  members  of  that  component. 

Lemma  19  Let  A  £  11  and  X  =  D(A)  —  A.  Let  be  the  current  matching. 
If  /i'  has  no  matches  between  members  of  X  and  members  of  A  and  is  not 
locally  perfect  on  A,  then  there  is  a  blocking  pair  (x,y)  for  p!  such  that  p,{x)  =  y. 

Proof.  Since  y'  is  not  a  locally  perfect  matching  on  A  there  must  be  some  xq  G  A 
such  that  y! (xq)  =  xg  or  yf  (xg)  G  S  —  X  —  A.  Let  yg  =  y{xg).  Now  since  y  is  a 
stable  matching,  yg  y'{xg).  If  xg  yy„  y'(yg)  then  (xo,?/o)  is  a  blocking  pair 
of  y'  and  y{xg)  =  yg. 

Otherwise  y' {yg)  Pyg  xg,  so  y'{yg)  G  D(^)-  In  fact,  y'  G  A  since  y'  has  no 
matches  between  members  of  A  and  X.  Let  xi  =  y' {yg)  and  yi  =  y{xi).  Since 
^  is  a  stable  matching,  yi  yo  or  else  {xi,yg)  would  form  a  blocking  pair  for 
y.  Now  if  xi  y'{yi),  {xi,yi)  is  a  blocking  pair  of  y'  and  y{xi)  =  yi  so  we 
have  our  result.  Otherwise  we  repeat  in  the  same  manner  to  form  a  sequence  of 
pairs  {{x„yi)}  such  that  Xi,yi  G  A,  y{xi)  =  y„  y' {yi)  =  Xj+i,  yi  y'{xf),  and 
Xi+i  Xi  for  all  i.  But  this  cannot  cycle  since  no  participant  is  repeated.  This 
is  because  at  each  step  we  add  a  new  pair  Xi,yi  where  y{xi)  =  yi  and  either 
y'{xg)  =  xg  or  y'{xg)  A,  so  Xg  cannot  be  repeated.  Furthermore,  it  cannot  go 
forever  since  A  is  finite.  Therefore  the  sequence  must  terminate  at  some  index  k 
and  {xk,yk)  is  a  blocking  pair  for  y' . 

Next  we  place  a  lower  bound  on  the  probability  that  we  make  some  progress 
toward  the  /i-stable  matching  when  a  stable  component  of  the  decomposition  is 
not  in  a  locally  perfect  matching. 

Lemma  20  Let  A  G  11  be  a  stable  component  of  size  at  most  c  and  X  =  D(A)  — 
A.  Let  y'  be  any  matching  on  S  that  is  not  a  locally  perfect  matching  on  A.  Then 
starting  from  y' ,  if  no  matches  are  formed  between  a  member  of  A  and  a  member 
of  X ,  the  probability  that  the  first  blocking  pair  resolved  between  two  members  of 
A  is  a  pair  in  y  is  at  least  ^ . 

Proof.  Lemma  19  shows  there  will  be  one  blocking  pair  which  is  in  y  until  the 
matching  becomes  locally  perfect  on  A.  In  order  for  the  matching  to  become 
locally  perfect  on  A,  a  blocking  pair  must  be  resolved  between  two  members  of 
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A.  Therefore  since  there  will  be  at  most  blocking  pairs  involving  two  members 
of  A  and  at  least  one  of  them  is  in  /x,  there  is  a  probability  that  the  first 
blocking  pair  resolved  between  members  of  A  is  in 

Using  this  lemma,  we  bound  the  probability  that  a  component  of  the  decom¬ 
position  will  make  some  progress  toward  the  /x-stable  matching  each  time  the 
matching  is  not  locally  perfect  on  it. 

Lemma  21  Let  A  €  11  be  a  stable  component  of  size  at  most  c  and  X  =  D(A)  — 
A.  Let  /io  be  any  matching  on  S  such  that  /xq  ]a  contains  m  of  the  pairs  in  /x 
where  0  <  m  <  c.  Let  yxo,/ii,  •  ■  •  ,/X(  be  any  valid  sequence  of  matchings  under 
the  better  response  dynamics  starting  from  /xq  such  that 

1.  Ht  is  locally  perfect  on  A 

2.  pLi  is  not  locally  perfect  on  A  for  some  i,0  <  i  <  t 

3.  /Xfe  does  not  have  any  matches  between  a  member  of  A  and  a  member  of  X 
for  some  k,0  <  k  <  t 

Then  the  probability  that  3j,  0  <  j  <  t,  qij  ]  a  contains  at  least  m  +  1  of  the  pairs 
in  pL  is  at  least 

Proof.  Assume  pto,  pti, . . . ,  pit  is  such  a  sequence,  and  i  is  the  first  index  such  that 
pii  is  not  locally  perfect.  Without  loss  of  generality  assume  k  =  t  is  the  hrst  index 
k  >  i  such  that  pi^  is  locally  perfect  on  A.  This  assumption  is  valid  because,  if 
there  is  at  least  a  probability  p  of  some  event  occurring  in  a  subsequence,  then 
there  is  clearly  at  least  a  probability  p  of  that  event  occurring  in  the  entire 
sequence. 

There  are  two  cases:  either  pig  is  locally  perfect  on  A  or  not. 
case  i:  Assume  /Xq  is  not  locally  perfect,  so  z  =  0.  Then  in  order  to  reach 
pit  there  must  be  at  least  one  match  formed  between  two  members  of  A.  Let 
j  >  0  be  the  hrst  index  in  the  sequence  such  that  pij  was  formed  by  resolving  a 
blocking  pair  between  two  members  of  A.  Since  no  one  in  A  prefers  anyone  in 
S'  — D(A)  to  his  partner  in  /x,  pij-i  ]a  has  m  pairs  in  pi.  By  lemma  20  there  is  at 
least  ^  probability  that  the  hrst  blocking  pair  resolved  between  two  members 
of  A  is  in  pi.  This  will  result  in  pij  ^a  having  m  J-  1  pairs  in  pi. 

case  ii:  If  pio  is  locally  perfect,  so  z  >  0.  There  are  two  ways  to  transition 
from  pii-i  to  pii.  One  is  for  a  blocking  pair  of  pii-i  between  a  member  of  A  and 
a  member  of  S  —  Ai  —  A  to  be  resolved.  Since  this  cannot  involve  a  member  of 
A  who  is  with  his  partner  in  pi  according  to  pi',  pit  ]a  has  m  pairs  that  are  in 
pi.  Therefore  this  case  reduces  to  the  hrst  case  where  the  initial  matching  is  not 
perfect. 

The  other  way  to  transition  from  /Xi_i  to  pit  is  for  a  blocking  pair  between 
two  members  of  A  to  be  resolved,  leaving  two  unmatched  members  of  A,  say  x,  y. 
The  blocking  pair  cannot  involve  two  pairs  of  pi  or  else  it  would  be  a  blocking 
pair  for  pi.  If  it  involves  no  pairs  of  pi  then  again  this  case  reduces  to  the  hrst 
case. 

In  the  last  case,  pit  ]  a  has  m  —  1  pairs  that  are  in  pi.  We  cannot  reach  pit 
without  resolving  a  blocking  pair  between  two  members  of  A.  Let  I  >  i  he  the 
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first  index  after  i  in  the  sequence  such  that  iii  was  formed  by  resolving  a  blocking 
pair  between  two  members  of  A.  Then  fii-i  ]  a  must  have  m  —  1  pairs  that  are  in 
fx.  By  lemma  20  there  is  at  least  ^  probability  that  /i;  has  m  pairs  that  are 
in  If  this  occurs,  the  blocking  pair  resolved  to  transition  to  /i;  cannot  involve 
both  X  and  y  because  they  are  not  partners  in  y,.  Thus,  at  least  one  of  a;  or  y  is 
still  not  matched  to  someone  in  A.  Therefore,  /x;  is  not  a  locally  perfect  matching 
on  A.  Then  by  the  first  case,  we  have  at  least  js  probability  that  for  some  j, 
I  <  j  <  t,  y,j  ]  A  has  TO  +  1  pairs  that  are  in  /x.  This  gives  us  a  total  probability 
of  at  least  ^  that  fj,j  ^a  has  to  +  1  pairs  that  are  in  y  for  some  j,  0  <  j  <t. 

We  now  bound  the  expected  number  of  times  each  stable  component  will 
have  to  become  not  locally  perfect  before  it  becomes  ^-stable. 

Lemma  22  Let  A  G  11  be  a  stable  component  of  size  at  most  c  and  X  =  D(A)  — 
A.  Let  p!  he  any  matching  on  S.  Then  starting  from  p' ,  if  no  matches  are  formed 
between  a  member  of  A  and  a  member  of  X,  the  expected  number  of  distinct 
times  the  matching  needs  to  transition  from  a  locally  perfect  matching  on  A  to 
a  matching  that  is  not  locally  perfect  on  A  before  it  reaches  a  p-stable  matching 
on  A  is  at  most  . 


Proof.  Consider  a  Markov  chain  with  states  {0, 1, ... ,  c}  where  state  i  represents 
a  matching  whose  restriction  to  A  has  i  pairs  in  p.  Let  U  be  the  expected  number 
of  times,  starting  from  state  i,  that  the  matching  transitions  from  a  locally 
perfect  matching  on  kl  to  a  matching  that  is  not  locally  perfect  on  A  before  it 
reaches  a  /x-stable  matching  on  A.  Then  tc  =  0  since  state  c  represents  a  /x-stable 
matching.  For  all  other  states,  by  lemma  21,  we  have  at  least  a  ^  probability 
of  reaching  state  i  +  1  from  state  i  after  one  or  fewer  transitions  from  a  locally 
perfect  matching  on  kl  to  a  matching  that  is  not  locally  perfect  on  A.  In  the 
worst  case,  we  will  move  to  state  0  after  one  such  transition  with  the  remaining 
probability.  This  leads  to  the  formula  ti  <  +  piti+i  +  1  for  x  =  0, 1, . . . ,  n. 

We  need  to  upper  bound  to  since  0  is  the  farthest  state  from  c.  Now  to  < 
ti  +  c'^.  Furthermore  if  to  <  U  +  ;  then 

_  1  ^  ^ 

to  <  — ^ — to  +  ^ti+i  +  1  +  ^  c'^^ 


so 


and 


■^to  <  -^U+i  +  ^  +  1 


i=i 


i+l 


to  ^  ^i+l  +  I  +  1  —  ti+l  +  I 


Aj 


G=1 


Therefore  to  <tc  +  =  J2j=i  ^ 


Aj  /  M(c+i) 
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The  final  lemma  we  need  shows  that  when  the  matching  is  not  locally  perfect 
on  a  stable  component  of  the  decomposition,  it  will  reach  a  perfect  matching 
in  time  that  depends  only  linearly  in  n  in  expectation,  provided  there  is  no 
interference  from  members  of  lower  stable  components. 

Lemma  23  Let  A  G  U  be  a  stable  component  of  size  at  most  c  and  X  =  D(A)  — 
A.  Let  p!  he  any  matching  on  S  which  is  not  locally  perfect  on  A.  Then  starting 
from  p' ,  if  no  matches  are  formed  between  a  member  of  A  and  a  member  of  X , 
the  expected  time  reach  a  matching  which  is  locally  perfect  on  A  is  at  most  cn^'^. 

Proof.  Lemma  19  implies  that  for  any  given  matching,  either  the  matching  is 
locally  perfect  on  A  or  there  is  a  blocking  pair  between  two  members  of  A  which 
is  a  pair  in  p.  Since  the  size  of  A  is  at  most  c,  there  are  at  most  c  such  pairs. 
Therefore  if  all  of  them  are  resolved  in  c  consecutive  steps,  the  resulting  matching 
will  be  locally  perfect  on  A.  Alternatively  if  after  fewer  than  c  steps  of  resolving 
blocking  pairs  that  are  in  p  we  reach  a  matching  with  no  such  blocking  pairs, 
then  the  matching  must  already  be  locally  perfect  on  A.  For  any  given  matching 
there  are  at  most  n'^  total  blocking  pairs  so  the  probability  of  resolving  a  blocking 
pair  between  two  members  of  A  that  is  a  pair  in  p  is  at  least  But  then  the 
probability  of  resolving  up  to  c  of  them  and  reaching  a  locally  perfect  matching 
in  c  or  fewer  steps  is  at  least 

Therefore,  in  expectation  we  will  have  to  repeat  the  process  of  making  c  steps 
at  most  times  before  reaching  a  locally  perfect  matching  on  A.  This  leads  to 
at  most  steps  in  expectation. 

Finally  we  will  show  that  the  expected  convergence  time  for  the  better  re¬ 
sponse  dynamics  is  linear  in  the  total  number  of  participants  but  possibly  ex¬ 
ponential  in  the  size  of  the  largest  stable  component  and  depth  of  the  decom¬ 
position.  The  special  case  where  the  size  of  the  decomposition  is  1  includes  the 
correlated  preferences  of  Ackermann  et  al. 

Theorem  24  (Convergence)  Suppose  p  is  a  stable  matching.  Suppose  the 
depth  of  {11,  is  d  and  the  size  of  the  largest  stable  component  of  11  is  no  more 
than  c.  Then  the  expected  time  to  converge  to  a  .stable  matching  is 
If  c=  1,  then  the  expected  time  is  0{n^). 

Proof.  Suppose  p'  is  another  stable  matching.  First,  suppose  that  for  any  stable 
component  A'  of  a  /i'-decomposition,  a  p'-stahle  matching  is  never  reached  on 

Consider  the  /^-decomposition  graph  for  (7T,  P).  Recall  that  a  stable  compo¬ 
nent  A  is  on  level  j  if  depth{A)  =  j.  For  convenience,  let  level  d -I- 1  be  an  empty 
dummy  level  at  the  top.  Since  the  depth  is  d,  there  are  exactly  d  -I-  1  levels. 
We  proceed  by  bounding  the  expected  time  for  one  level  to  reach  a  p  stable 
matching,  and  then  recurse  on  the  higher  levels. 

Let  T{1)  denote  the  expected  time  for  the  participants  in  stable  components 
on  levels  I  and  above  to  reach  a  stable  matching  without  resolving  blocking 
pairs  involving  any  members  of  stable  components  on  lower  levels.  Let  ni  be  the 
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number  of  stable  components  on  level  1.  Note  that  since  there  are  at  most  n  stable 
components  of  U,  ni  + . .  .  +  nd  <  n.  We  will  show  that  T{0)  = 

First  observe  that  T{d+  1)  =  0  since  there  are  no  stable  components  at  level 
c?  +  1. 

Now  consider  T{1)  for  ?  <  c?  +  1. 

When  one  of  the  n/  stable  components  A  on  level  I  is  not  in  a  locally  perfect 
matching.  Then  by  Lemma  23,  we  know  it  will  take  steps  in  expectation 
to  reach  a  locally  perfect  matching  on  A.  Also,  by  lemma  21  we  know  it  has  at 
least  ^  probability  of  reaching  a  matching  whose  restriction  to  A  has  a  greater 
number  of  pairs  that  are  in  /r  than  the  current  matching,  before  it  reaches  a 
locally  perfect  matching. 

On  the  other  hand,  when  all  ni  stable  components  are  in  locally  perfect 
matchings,  then  there  are  two  cases: 

If  there  is  a  blocking  pair  between  two  members  of  stable  components  on 
level  I  it  will  remain  there  until  the  matching  becomes  not  locally  perfect  on 
at  least  one  stable  component  on  level  1.  Since  there  are  at  most  blocking 
pairs,  it  will  take  at  most  steps  in  expectation  for  the  matching  to  become 
not  locally  perfect  on  at  least  one  stable  component  on  level  1. 

If  there  are  no  such  blocking  pairs,  it  might  be  required  for  the  higher  levels  to 
reach  a  stable  matching  before  exposing  a  blocking  pair  involving  a  participant 
on  level  1.  If  no  matches  are  formed  involving  any  members  of  components  on 
level  I  or  lower,  the  expected  time  for  the  remaining  stable  components  to  reach 
a  stable  matching  is  given  by  T{1  +  1).  Once  the  higher  levels  have  reached  a 
stable  matching,  the  only  blocking  pairs  not  involving  members  of  levels  below  I 
are  between  a  member  of  a  stable  component  on  level  I  and  a  member  of  a  stable 
component  on  a  higher  level.  Unless  all  stable  components  on  level  I  and  above 
are  in  a  stable  matching,  at  least  one  such  blocking  pair  must  exist.  Therefore 
it  will  only  take  1  more  step  to  reach  a  matching  which  is  not  locally  perfect  on 
one  stable  component  on  level  1. 

Consequently,  it  will  take  at  most  +  r(^  +  1)  +  1  steps  to  reach  a  matching 
that  is  not  locally  perfect  on  one  stable  component  on  level  1.  Again,  by  Lemma 
23,  we  know  it  will  take  cn^'^  steps  in  expectation  to  reach  a  locally  perfect 
matching  on  A.  By  Lemma  22,  we  know  in  expectation,  for  each  stable  com¬ 
ponent  on  level  I,  it  will  take  at  most  c'^C+i)  transitions  from  a  locally  perfect 
matching  to  a  matching  which  is  not  locally  perfect  on  that  stable  component  it 
reaches  a  /r-stable  matching.  This  means  that  in  expectation  it  will  take  at  most 
of  these  transitions  total  before  all  stable  components  on  level  I  reach 
a  /x-stable  matching. 

Therefore,  in  the  worst  case,  it  will  take  (n^  -|-T(/ -b  1)  -I- 1)  steps  to  transition 
from  a  locally  perfect  matching  to  a  matching  that  is  not  locally  perfect  on  one 
of  the  stable  components  on  level  1.  Then  it  will  take  at  most  cn?‘^  steps  to  reach 
a  matching  which  is  locally  perfect  on  that  stable  component.  Furthermore,  this 
process  needs  to  be  repeated  no  more  than  times  in  expectation  in 

order  for  all  stable  components  on  level  I  to  reach  a  /i-stable  matching. 
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Once  all  stable  components  on  level  I  have  reached  a  /x-stable  matching,  all 
that  remains  is  for  the  higher  levels  to  reach  a  stable  matching,  which  takes 
T{1  +  1)  time  in  expectation. 

This  yields  the  following  formula: 


Til)  <  niC^(“+i)(cn2"  +  n2+r(Z  +  l)  +  l)+T(/  +  l)  <  +  Til +  1)) 

Solving  this  recursion  for  T(0),  we  obtain 

r(0)  <  2noc‘‘(‘=+i)(cn2=  +  T(1)) 

d+1  i-1 

i=l  j=0 

so  since  +  1  <  0(n)  for  all  i,  r(0)  = 

This  is  the  expected  time  to  reach  the  stable  matching  /i.  Now  suppose  for 
some  stable  component  of  a  /i'-decomposition  for  some  other  stable  matching 
/i',  a  /r'-stable  matching  is  reached  on  D^/(A').  By  Lemma  18,  this  will  not 
increase  the  size  or  depth  of  the  remaining  decomposition.  Therefore,  if  this 
happens  before  /r  is  reached,  it  will  only  decrease  the  convergence  time. 

Finally,  as  a  special  case  assume  c  =  1.  In  this  case  a  locally  perfect  matching 
on  a  stable  component  is  a  /i-stable  matching.  By  lemma  23  it  will  take  at  most 
steps  for  a  stable  component  on  level  I  to  reach  a  /x-stable  matching.  Since 
there  are  n;  components  on  level  /,  T{1)  <  m'n?  +  T{1  +  1)  < 

^(0)  =<  Y.i=i  =  n^- 

4  Correlated  and  Intercorrelated  Preferences 

We  have  shown  bounds  on  convergence  time  but  this  is  only  relevant  if  there 
is  variation  in  the  jealousy  graph  structures  of  real  markets.  While  randomly 
generated  preferences  tend  to  have  decompositions  that  are  close  to  the  trivial 
decomposition,  which  is  the  entire  set,  real-world  markets  tend  to  have  some 
structure.  Here  we  show  that  two  classes  of  preferences  found  in  real  world 
markets,  correlated  and  intercorrelated  preferences,  exhibit  decompositions  with 
small  size  components.  Partially  correlated  preferences  are  often  used  by  model¬ 
ers  [8, 11]  and  are  natural  in  many  matching  markets  (e.g.  mate  selection)  where 
preferences  are  based  on  a  mixture  of  universally  desirable  features  (e.g.  intel¬ 
ligence)  and  idiosyncratic  tastes  (e.g.  shared  hobbies).  Note  that  the  correlated 
preferences  discussed  here  differ  from  the  correlated  preferences  of  Ackemann 
et  al.  Intercorrelation  exists  when  the  preferences  of  the  men  relate  to  the  pref¬ 
erences  of  the  women.  See  [12]  for  examples  of  markets  with  intercorrelation. 
Boudreau  showed  that  more  correlation  and  intercorrelation  lead  to  faster  con¬ 
vergence  of  the  better  response  algorithm  [8].  We  provide  similar  plots  in  Figures 
1(b)  and  1(d).  Theorem  24  provides  theoretical  justihcation  for  these  simulated 
results. 


Jealousy  Graphs 


13 


Correlation  Factor  (U) 


(a)  Correlation  vs.  Size  and  Depth 


(b)  Correlation  vs.  Convergence  Time 


Intercorrelation  Factor  (V) 
(c)  Intercorrelation  vs.  Size 
Depth 


Intercorrelation  Factor  (V) 


(d)  Intercorrelation  vs.  Convergence 
Time 


Fig.  1.  Jealousy  Graphs  vs.  Correlation  and  Intercorrelation,  (a)  The  jealousy  graph 
parameters  change  as  preferences  become  more  correlated,  (b)  Convergence  time  de¬ 
creases  as  preferences  become  more  correlated,  (c)  The  jealousy  graph  parameters 
change  as  preferences  become  more  intercorrelated.  (d)  Convergence  time  decreases  as 
preferences  become  more  intercorrelated. 


As  described  in  [13, 11],  correlated  preferences  are  generated  using  scores  of 
the  form: 


Synw  —  Vmw  U 


where  S^w  is  the  score  man  m  gives  woman  w  composed  of  his  individual  score 
Vmw  and  a  correlation  factor  U  €  [0,oo)  multiplied  by  the  consensus  score  of 
w,  Iw  ilmw  and  ly,  are  chosen  uniformly  at  random  from  [0,  Ij.  The  men  then 
rank  the  women  in  order  from  lowest  score  to  highest.  Women’s  preferences  are 
generated  analogously.  For  various  values  of  U  we  generate  100  preferences  with 
correlation  factor  U.  For  each  set  of  preferences  we  find  the  decomposition  with 
smallest  size  and  report  the  average  of  these  sizes.  We  also  compute  the  average 
minimal  depth  in  the  same  manner.  The  results  are  shown  in  figure  1(a).  At 
U  =  0,  the  average  size  is  close  to  n  and  the  depth  is  close  to  1.  As  17  goes  to 
oo,  the  average  size  approaches  1  and  the  depth  approaches  n.  These  are  the 
parameters  of  perfectly  correlated  preferences.  This  shows  that  as  the  amount  of 
correlation  varies,  so  do  the  size  and  depth  of  the  decompositions.  Figure  1(b) 
shows  the  log  of  the  average  convergence  time  over  100  trials  for  each  of  the  100 
correlated  preferences  generated. 
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As  in  [12],  intercorrelated  preferences  can  be  generated  using  scores  of  the 
form: 


SrriiWj  —  VmiWj  “t“  *  ji  j\n 

where  Srmwj  is  the  score  man  rrii  gives  woman  Wj.  As  with  correlated  pref¬ 
erences,  rjrrnwj  is  his  individual  score.  Here  V  is  the  intercorrelation  factor  and 
I*  ~  j\n  =  min(|i  —  j\,n  —  |*  —  j|)/(f )  represents  the  “distance”  man  rrii  is  from 
woman  Wj.  1(c)  and  l((i)  are  generated  in  the  same  manner  as  1(a)  and  1(&), 
respectively.  These  plots  show  that  as  preferences  become  more  intercorrelated, 
the  size  and  depth  of  the  decompositions  decrease.  As  Theorem  24  explains,  this 
decreases  the  convergence  time  of  the  better  response  algorithm  as  intercorrela¬ 
tion  increases. 

5  Conclusion 

We  have  introduced  a  new  way  of  viewing  stable  matching  problems  in  terms  of 
their  jealousy  graphs  and  /^-decompositions.  We  demonstrate  that  these  concepts 
are  useful  in  analyzing  the  convergence  time  of  the  better  response  algorithm 
and  guarantee  polynomial  convergence  on  a  subclass  of  matching  markets.  Fur¬ 
thermore,  these  theoretical  results  apply  to  a  broad  range  of  markets  since  they 
provide  a  notion  of  structure  which  extends  beyond  the  well-studied  notions  of 
correlation  and  intercorrelation. 

One  open  question  involves  the  exponential  dependency  on  the  depth  of  the 
decomposition.  While  we  know  that  the  exponential  dependency  on  size  cannot 
be  removed,  it  remains  an  open  question  whether  we  can  improve  this  bound  in 
terms  of  the  depth.  Another  open  problem  concerns  which  matching  is  most  likely 
to  be  reached.  Since  our  result  provides  a  method  of  classifying  the  expected  con¬ 
vergence  time  of  the  better  response  algorithms  in  terms  of  the  decompositions 
of  the  stable  matchings,  we  conjecture  that  matchings  with  decompositions  that 
have  small  size  and  depth  are  more  likely  to  be  reached  than  ones  with  large  size 
and  depth. 
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Appendix 

Proof  of  Theorem  14 

Proof.  Let  G  =  {V,E)  and  S  =  {mi,Wi\i  =  1,2, ...  ,n}.  Arbitrarily  index  the 
vertices  in  U  as  Ui,  V2,  ■  •  ■ ,  v„.  We  will  design  our  preferences  such  that  if  = 

Wi  for  f  =  1,  2, . . . ,  n,  /i  is  a  stable  matching  and  Vi  is  the  vertex  corresponding 
to  {mi,Wi}.  Now  define  P  =  {PmuPwi  |1  <  z  <  n}  as  follows.  First,  for  every 
woman  Wi,  define  such  that  rrii  rrij  for  j  ^  i.  The  remaining  ordering 
can  be  arbitrary.  For  every  man,  define  such  that  wj  Prai  Wi  44  (vj,Vi)  G  E. 
The  ordering  among  elements  within  {wj\wj  y-mi  Wi}  and  {wj\wi  ymi  Wj}  can 
be  arbitrary. 

To  see  that  ^  is  indeed  a  stable  matching,  observe  that  under  /r  all  women  are 
matched  with  their  top  choice.  Therefore,  no  woman  has  incentive  to  deviate,  so 
there  can  be  no  blocking  pairs  and  /x  is  a  stable  matching  on  S,  P. 

All  that  remains  is  to  show  =  G.  Now  in  =  (V',E')  let  the  vertices 
be  denoted  V'  =  . . . ,  We  will  let  v'  correspond  to  {mi,Wi}  for  i  = 

1,2, ...  ,n.  Since  all  women  are  matched  with  their  top  preference  by  /x,  the 
women  will  not  be  responsible  for  any  edges  in  J^.  Therefore,  (u',z;')  G  A'  <f4 
Wi  yrrij  Wj  <t4  (vi,Vj)  G  A  by  the  way  we  defined  for  all  i,j.  Thus  if  we  let 
Vi  =  v[  we  have  equivalent  graphs. 


Proof  of  Lemma  15 

Proof.  Let  vi,V2, .  ■ .  ,Vk  be  the  vertices  along  the  path  such  that  vi  = 
and  Vk  =  v^m'.w'}-  Let  Ai  be  the  stable  component  containing  Vi.  Then  since  the 
partial  order  is  induced  by  the  edges  between  components  of  J^,  either  Ai  = 
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or  Ai  :<  Ai+i.  Therefore,  by  transitivity  Ai  <  Ak  =  A.  Thus  m,w  G  D(A) 
because  V{rn,w}  G  Gli. 

Proof  of  Lemma  16 

Proof.  Let  A  G  11  and  a  G  A.  Suppose  there  is  some  s  G  S—D(A)  such  that  s 

Then  since  s  ^  A,  there  are  distinct  vertices  in  J^,  Va,Vs  corresponding  to 
the  pair  with  a  and  the  pair  with  s,  respectively.  Edge  {vs,Va)  must  also  be  in 
since  s  )^a  /r(a).  Thus  there  is  a  path  in  from  Vg  to  Va,  so  by  Lemma  15, 
sGT>{A).  This  is  a  contradiction. 


Proof  of  Lemma  17 

Proof.  Suppose  there  is  x  G  but  fj,{x)  G  A  —  B.  By  Lemma  16,  f.L'{x) 

fj.{x).  In  that  case,  fJ,'{x)  G  Ap|iJ,  also  by  Lemma  16.  Now  since  x  prefers 
/i'(x)  to  /i(x),  /x(/i'(x))  X  or  else  {x,^'{x))  forms  a  blocking  pair  for 

Again  by  Lemma  16,  fj.{fj.'{x))  G  A(^B.  Continuing  in  this  manner  gives  an 
infinite  sequence  x ,  fi' (x) ,  (x)) ,  fj,' (x))) , . . .  G  Ap|i3.  These  elements  are 

distinct  elements  since  ^{x)  fJ-'{x)  and  both  p.  and  fj,'  are  bijective.  This  is  a 

contradiction  since  Ap|i?  is  finite.  Therefore  /x(Ap|iJ)  =  fj,'{Af]B)  =  A(^B. 


Proof  of  Lemma  18 

Proof.  Let  be  such  that  /x"  1d(A)=  1d(A)  and  /x"  1s_d(.4)=  M  1s-d(A)- 

Clearly  there  are  no  blocking  pairs  involving  two  members  of  D(A)  or  else  /i' 
would  not  be  stable  and  there  are  no  blocking  pairs  between  two  members  of 
S  —  D(A)  or  else  /x  would  not  be  stable.  Finally,  by  Lemma  16  no  member  of 
D(A)  can  prefer  any  member  of  S'  —  D(A)  to  his  partner  in  /x'.  Therefore  there 
can  be  no  blocking  pairs  between  a  member  of  D(A)  and  a  member  of  S  —  D(A) 
so  is  indeed  a  stable  matching. 

By  Lemma  17,  for  each  stable  component  B  of  /x(i3  — D(A))  =  B  —  T>{A). 
The  set  {B  —  D(A)|i3  G  p^}  forms  a  partition  of  S  —  D(A)  and,  paired  with 
the  same  partial  order  as  p^,  forms  a  decomposition.  Clearly  the  size  has  not 
increased  since  the  sets  in  the  partition  are  no  larger  and  the  depth  has  not 
increased  since  the  decomposition  has  the  same  partial  order  as  p^. 


